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Abstract 
Background: Human infections with zoonotic coronaviruses (CoVs), including severe acute respiratory syndrome (SARS)-CoV and 
Middle East respiratory syndrome (MERS)-CoV, have raised great public health concern globally. Here, we report a novel bat- 
origin CoV causing severe and fatal pneumonia in humans. 

Methods: We collected clinical data and bronchoalveolar lavage (BAL) specimens from five patients with severe pneumonia from Jin 
Yin-tan Hospital of Wuhan, Hubei province, China. Nucleic acids of the BAL were extracted and subjected to next-generation 
sequencing. Virus isolation was carried out, and maximum-likelihood phylogenetic trees were constructed. 

Results: Five patients hospitalized from December 18 to December 29, 2019 presented with fever, cough, and dyspnea accompanied 
by complications of acute respiratory distress syndrome. Chest radiography revealed diffuse opacities and consolidation. One of 
these patients died. Sequence results revealed the presence of a previously unknown B-CoV strain in all five patients, with 99.8% to 
99.9% nucleotide identities among the isolates. These isolates showed 79.0% nucleotide identity with the sequence of SARS-CoV 
(GenBank NC_004718) and 51.8% identity with the sequence of MERS-CoV (GenBank NC_019843). The virus is phylogenetically 
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closest to a bat SARS-like CoV (SL-ZC45, GenBank MG772933) with 87.6% to 87.7% nucleotide identity, but is in a separat 
clade. Moreover, these viruses have a single intact open reading frame gene 8, as a further indicator of bat-origin CoVs. However, the 
amino acid sequence of the tentative receptor-binding domain resembles that of SARS-CoV, indicating that these viruses might use 


the same receptor. 


Conclusion: A novel bat-borne CoV was identified that is associated with severe and fatal respiratory disease in humans. 
Keywords: Bat-origin; Coronavirus; Zoonotic transmission; Pneumonia; Etiology; Next-generation sequencing 





Introduction 


Coronaviruses (CoVs) are enveloped viruses with a single 
positive-stranded RNA genome (~26-32 kb in length). 
They belong to the subfamily Orthocoronavirinae under 
the family Coronaviridae, and are classified into four 
genera: Alphacoronaviruses (a), Betacoronaviruses (B) 
Gammacoronaviruses (y), and Deltacoronaviruses (8).'” | 
The viral genome normally encodes four structural 
proteins, spike (S), envelope (E), membrane (M), and 
nucleocapsid (N), as well as several non-structural proteins 
and multiple unique accessory proteins.!!! 


CoVs infect humans and a variety of avian and mammalian 
species worldwide. There are six CoVs known to infect 
humans, including two a-CoVs (229E and NL63) and four 
B-CoVs (OC43, HKU1, severe acute respiratory syndrome 
[SARS]-CoV, and Middle East respiratory syndrome 
[MERS]-CoV).!!4) All human CoVs are zoonotic as a 
distinguishing characteristic.! In particular, bats are 
regarded as a key reservoir of CoVs, and many human 
CoVs are believed to have originated from bats.'°*! Since 
the beginning of this century, two zoonotic CoVs, SARS- 
CoV and MERS-CoV, have been identified to cause severe 
human diseases.'*-*7! The outbreak of SARS-CoV in 2003 
was responsible for 8096 cases and 794 deaths world- 
wide.'”' Since its discovery in Middle Eastern countries in 
2012, MERS-CoV has infected 2260 people with a current 
case fatality rate of 35.5%.'?1°! These outbreaks have 
raised public health concerns of the potential for the 
emergence of another novel zoonotic CoV. 


Here, we report a previously unknown bat-origin CoV 
causing severe and fatal pneumonia in five patients from 
Wuhan, China. Sequence results revealed that this virus, 
harboring a single open reading frame gene 8 (ORF$8), is 
phylogenetically closest to bat SARS-like CoV, but is in a 
separate lineage. Furthermore, the amino acid sequence of 
the tentative receptor-binding domain (RBD) of this new 
CoV resembles that of SARS-CoV, indicating that they 
might use the same receptor. These findings highlight the 
urgent need for regular surveillance of the interspecies 
transmission of bat-origin CoV to human populations. 


Methods 


Ethical approval 


This study was conducted in accordance with the 
Declaration of Helsinki and was approved by the National 
Health Commission of China and Ethics Commission of 
the Jin Yin-tan Hospital of Wuhan (No. KY-2020-01.01). 
The requirement for written informed consent was waived 
given the context of emerging infectious diseases. 


Clinical specimen and data collection 


Bronchoalveolar lavage fluid (BAL) samples were collected 
from five patients hospitalized with pneumonia in Jin Yin- 
tan Hospital of Wuhan, Wuhan, Hubei province, China 
from December 18 to 29, 2019. Information was gathered, 
including clinical data, demographic characteristics, 
underlying medical conditions, clinical signs and symp- 
toms, chest radiographic findings, clinical laboratory 
testing results, traveling history, recent animal exposure, 
and outcomes. The data collected for the cases were 
deemed by the National Health Commission of the 
People’s Republic of China as the contents of a public 
health outbreak investigation. 


Genome sequencing 


Nucleic acids were extracted from 200 uL BAL of each 
sample with the Direct-zol RNA Miniprep kit (Zymo 
Research, Irvine, CA, USA) and Trizol LS (Thermo Fisher 
Scientific, Carlsbad, CA, USA) according to the manu- 
facturer instructions’s in a biosafety III laboratory. A 50- 
wL elution was obtained from each sample. The DNA/ 
RNA concentrations were measured by a Qubit Fluo- 
rometer (Thermo Fisher Scientific). The sequencing 
library was constructed by a transposase-based method- 
ology and sequenced on an Illumina sequencing platform 
(Illumina, San Diego, CA, USA). At least 25 million 
single-end 76-bp reads were generated for each sample on 
the Illumina NextSeq platform. Quality control processes 
included removal of low-complexity reads by bbduk 
(entropy =0.7, entropy-window =50, entropy k =5; 
version: January 25, 2018),11) adapter trimming, low- 
quality reads removal, short reads removal by Trimmo- 
matic (adapter: TruSeq3-SE.fa:2:30:6, LEADING: 3, 
TRAILING: 3, SLIDING WINDOW: 4:10, MINLEN: 
70, version: 0.36),''7! host removal by bmtagger (using 
human genome GRCh38 and yh-specific sequences as 
reference),'1?! and ribosomal reads removal by Sort- 
MeRNA (version: 2.1b).''4! Taxonomic assignment of 
the clean reads was performed with Kraken 2 against the 
reference databases, including archaea, bacteria, fungi, 
human, plasmid, protozoa, univec, and virus sequences 
(software 2.0.7-beta, database version: August 2, 
2019).!'5! A negative control sample was processed and 
sequenced in parallel for each sequencing run as a 
contamination control. The data were classified by 
simultaneous alignment to the microbial genome data- 
bases comprising viruses, bacteria, fungi, and parasites 
after filtering of the adapters and human-origin reads. 
The sequences were confirmed by Sanger sequencing 
with specific primers and one-step real-time polymerase 
chain reaction (RT-PCR) Kit (Invitrogen, Carlsbad, 
CA, USA). 
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Phylogenetic analysis 


Multiple sequence alignment was performed with the 
ClustalW program using MEGA software (version 
7.0.14). Phylogenetic trees were constructed by means of 
the maximum-likelihood method with MEGA software 
(version 7.0.14). The full-genome viral sequences were 
deposited in the dataset of Global Initiative on Sharing All 
Influenza Data (No. EPIISL_402123, EPI_ISL_403928- 
31) and the Genome Warehouse in National Genomics Data 
Center, Beijing Institute of Genomics, Chinese Academy of 
Sciences, under Project ID PRJCA002165 that is publicly 
accessible at https://bigd.big.ac.cn/gwh as of January 2020. 


Virus isolation 


The BAL specimens were inoculated onto Vero cells 
(American Tissue Culture Colection [ATCC], CCL-81). 
All cultures were observed daily for a cytopathic effect 
(CPE). Maintenance medium containing tosyl-phenylala- 
nine chloromethyl-ketone enzyme ata final concentration of 
1 pg/mL was replenished at day 4, and cultures were 
terminated 7 days after inoculation. The viral particles were 
negative stained with 1% solution of phosphotungstic acid 
(pH 7.0) and the morphology was characterized by using 
120 kV TECNAI (Thermo Fisher Scientific, Hillsboro, OR, 
USA) electron microscopy and camera of Gatan832 (Gatan, 
Pleasanton, CA, United States). The culture supernatants of 
cells demonstrating CPE were mixed with paraformalde- 
hyde, dried onto formvar/carbon-coated grids, and stained. 
Viral nucleic acids were confirmed by RT-PCR with specific 
primers [Supplementary Table 1, http://links.lww.com/ 
CM9/A191]. 


Immunofluorescence assay 


Spot slides were prepared by applying 20 wL of the virus- 
infected or non-infected cell suspension onto 12-well 
Teflon-coated slides. The cells were fixed with 4% 
paraformaldehyde in 1x  phosphate-buffered saline 
(PBS) for 30 min, washed three times with PBS, blocked, 
and stained with serum from a convalescent patient or 
serum from a healthy person for 30min at 37°C at a 
dilution of 1:200. Goat anti-human immunoglobulin G 
conjugated with fluorescein isothiocyanate was used as the 
secondary antibody (Jackson Immuno Research Labora- 
tories, Inc., West Grove, PA, USA). Nuclei and the 
cytoplasm were counterstained with 4’,6’-diamidino-2- 
phenylindole and Evans blue (Sigma-Aldrich, St. Louis, 
MO, USA). Fluorescent images were obtained and 
analyzed using laser-scanning confocal microscopy (Air- 
yscan LSM880, Zeiss, Berlin, Germany). 


Results 


General information of patients 


Patient 1 was a 65-year-old man who reported a high fever 
and cough, with little sputum production, at the onset of 
illness. He had a continuous fever and developed severe 
shortness of breath 16 days later. He was a vendor at the 
Huanan Seafood Market, Wuhan, Hubei Province, China. 
Patient 2, a 49-year-old woman, presented with high fever 
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and dry cough. Five days later, she developed dyspnea and 
was admitted to the hospital. She was also a worker in the 
Huanan Seafood Market. Patient 3 was a 52-year-old 
woman who did not report any market exposure. She was 
admitted to hospital because of fever, cough, and ground- 
glass opacity in the chest computed tomography scan. 
Patient 4 was a 41-year-old man who also presented with 
high fever and dry cough at the onset of the illness. He 
developed acute respiratory distress syndrome 7 days later. 
This patient had no known history of exposure to the 
Huanan Seafood Market. Patient 5, a 61-year-old man, 
was admitted to a local hospital with a 7-day history of 
fever, cough, and dyspnea. He also worked in the market. 


With regards to medical history, Patient 4 had hyperten- 
sion, and Patient 5 had chronic liver disease and abdominal 
myxoma, whereas none of the other patients had a record 
of underlying diseases. The demographic and clinical 
characteristics of the five patients are summarized in 
Table 1. 


Novel CoV identification by next-generation sequencing 


The resultant clean reads accounted for 12.0% to 92.0% 
of the raw reads. Most of the reads could be successfully 
assigned. Notably, 80.3% of the reads mapped to the viral 
genome for sample from Patient 5 with the highest 
proportion of viral reads among the five samples. Nearly 
all of the viral reads (97%) were classified as Coronavir- 
idae. Similarly, in the other four patients, most of the viral 
reads were assigned to B-CoVs. Based on de novo assembly 
and careful curation, a consensus sequence of this CoV was 
obtained. 


A substantial proportion of all sequencing reads mapped to 
the newly reported CoV genome (BWA mem, version: 
0.7.12), "6l ranging from 71,883 (0.27% among all reads) 
in Patient 4 to 37,247,818 (85%) in Patient 5. In addition, 
very few reads mapped to known bacterial pathogens, 
including Streptococcus, Acinetobacter baumannii, and 
Pseudomonas [Figure 1A-F]. 


The reads mapping to CoVs were assembled, and their 
genome sequences were confirmed by Sanger sequencing. 
The nucleotide (nt) similarity among the obtained five 
whole-genome sequences was 99.8% to 99.9%. The full 
length of the obtained genome was 29,870 bp with a GC 
content of 37.99% to 38.02%. The genome organization, 
5'-ORF1lab-S-E—M-N-3’, was similar to that of the most 
well-known bat SARS-like (SL)-CoV [Figure 2A]. In 
addition, unique accessory ORFs were identified that are 
characterized in the subgenus Sarbecovirus, encoding 
putative ORF3, ORF6, ORF7, and ORFS proteins reading 
from the 5’-terminus to the 3’-terminus between the 
structural proteins [Figure 2A]. 


Phylogenetic analysis 


Homology assessment showed that full-length viral genome 
sequences have 79.0% nt identity with that of SARS-CoV 
Tor2 (GenBank NC_004718), 51.8% with that of MERS- 
CoV (GenBank NC_019843), and 87.6% to 87.7% with 
those of bat SL-CoV ZC45 and ZXC21 (GenBank 
MG772933, MG772934), isolated from Chinese horseshoe 
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Table 1: Demographic, epidemiologic, clinical presentations, and treatment of the patients. 





Characteristics Patient 1 Patient 2 Patient 3 Patient 4 Patient 5 

Age, year 65 49 52 41 61 

Gender Male Female Female Male Male 
Huanan Seafood Market exposure Yes Yes No No Yes 
Smoking No No No No No 
Underlying conditions No No No Hypertension Chronic liver disease, 


Date of illness onset 
Date of admission 
Admission to ICU 
Sampling date 
Signs and symptoms 
Cough 
Sputum production 
Dyspnea (onset date) 
Fever 
Muscle pain/fatigue 
Headache 
Diarrhea 
Highest temperature, °C 
Laboratory tests on admission 
Number of white blood cells, x 10°/L 
Number of neutrophils, x 107/L 
Number of lymphocytes, x 10°/L 
Number of platelets, x 10°/L 
Prothrombin time, s 
APTT, s 
D-dimer, mg/L 
ALT, U/L 
Serum potassium, mmol/L 
Serum creatine, pmol/L 
Procalcitonin, ng/mL 
Pathogen testing 
Coronavirus reads proportion, % 
Imaging 
Bilateral GGO 
Consolidation 
Treatment 
Oxygen therapy 


ECMO (initiation date) 
Antibiotic therapy 
CRRT 
Complications 
ARDS (onset date) 
Septic shock (onset date) 
AKI (onset date) 
Secondary infections 
Outcome’ 


Dec 15, 2019 
Dec 18, 2019 
Dec 22, 2019 
Dec 24, 2019 


Yes 
Dec 31, 2019 


11.9 
11.6 
0.2 
92 
12 
25.5 
40.7 
50 
5.0 
53.1 
1.46 


39.9 


Yes 
Yes 


Mechanical 
ventilation 
Jan 6, 2020 
Yes 
No 


Dec 22, 2019 


Hospitalized 


Dec 22, 2019 
Dec 27, 2019 
Dec 29, 2019 
Dec 30, 2019 


Yes 
No 

Dec 27, 2019 
Yes 


0.44 
273 
12.5 
25.2 
0.21 
45 
3.8 
42.6 
<0.05 


73.1 


Yes 
Yes 


HFNC 


Dec 29, 2019 


Hospitalized 


Dec 22, 2019 
Dec 29, 2019 
Dec 29, 2019 
Dec 30, 2019 


Yes 
No 
Not available 


140 
12.3 
28.5 
0.65 
19 
4.0 
54.5 
<0.05 


13.6 


Yes 
No 


Nasal 
cannula 
No 
Yes 


No 
Discharged 


Dec 16, 2019 
Dec 22, 2019 
Dec 30, 2019 
Dec 30, 2019 


Yes 
No 
Dec 22, 2019 
Yes 
Yes 
No 
No 
39 


6.6 
5.0 
0.98 
129 
11.2 
26.9 
1.31 
29 
5.4 
85.1 
<0.05 


1.6 


Yes 
Yes 


Non-invasive 
ventilation 


Hospitalized 


myxoma of 
abdominal cavity 
Dec 20, 2019 
Dec 27, 2019 
Dec 28, 2019 
Jan 1, 2020 


Yes 
No 
Dec 27, 2019 

Yes 

Yes 

No 

No 
N.A. 


17.9 
16.2 
1.2 
315 
28.4 
105.2 
20.6 
35 
4.8 
106.6 
0.69 


80.3 


Yes 
Yes 


Mechanical 
ventilation 
Jan 2, 2020 
Yes 
Yes 


Dec 27, 2019 
Dec 31, 2019 
Dec 31, 2019 
Yes 
Death 


* Outcome data as of January 11, 2020. ICU: Intensive care unit; APTT: Activated partial thromboplastin time; ALT: Alanine aminotransferase; GGO: 
Ground-glass opacity; HFNC: High-flow nasal cannula; ECMO: Extracorporeal membrane oxygenation; CRRT: Continuous renal replacement 
therapy; ARDS: Acute respiratory distress syndrome; AKI: Acute kidney injury. 


bats (Rhinolophus sinicus) [Table 2], indicating that the 


novel CoVs are most similar to bat SL-CoVs. 


Compared with bat SL-CoV ZC45, the novel CoVs 
showed 75.9%, 98.6%, 93.2% to 93.4%, and 91.1% nt 


identities in the S, E, M, and N genes, respectively. Overall, 


ORF1lab showed 89.0% nt identity between the novel 
CoVs and bat SL-CoV ZC45. Surprisingly, RNA-depen- 


dent RNA polymerase (RdRp), which is the most highly 
conserved sequence among different CoVs,!'*! 


only 
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= Betacoronavirus 

= Acinetobacter 

= Pseudomonas 

= Serratia 

= Delftia 

= Yarrowia 

= Rhodococcus 

= Brevundimonas 

= Methylobacterium 

= Stenotrophomonas 

= Cutibacterium 

= Bradyrhizobium 

= Sphingomonas 
Others 


= Escherichia 

= Streptococcus 

= Betacoronavirus 

= Gemella 

= Lactobacillus 
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= Salmonella 
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= Methylobacterium 
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= Betacoronavirus 

= Chryseobacterium 
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= Enterococcus 
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= Lactobacillus 
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= Mycoplasma 
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a Brevibacterium 
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Others 





Figure 1: Microbial species profiles (proportions) in bronchoalveolar lavage fluid specimens from Patient 1 (A), Patient 2 (B), Patient 3 (C), Patient 4 (D), and Patient 5 (E) analyzed by deep 


sequencing. 


showed 86.3% to 86.5% nt identities with bat SL-CoV 
ZC45. According to the International Committee on 
Taxonomy of Viruses criteria, a new CoV species could be 
defined if the nt identity is less than 90% for the conserved 
RdRp sequence.'*! Thus, we considered that the novel 
CoVs should be classified as a new species under the 
subgenus Sarbecovirus of the genus Betacoronavirus. 


The phylogenetic trees constructed with the sequences of 
the RdRp, S, and N genes, and the whole genome using a 
maximum-likelihood model showed that all five novel 
CoVs were closely related to bat SL-CoVs ZXC21 and 
ZC45, but in a separate evolutionary lineage under the 


subgenus Sarbecovirus [Figure 2B-E], which is consistent 
with the homology assessment results. 


ORF3 and intact ORF8 gene regions were present in the 
novel CoVs, which are the characteristic features of bat- 
origin CoVs.!'7'8! ORF3 of the novel CoVs showed 
87.8% nt and 90.9% amino acid (aa) identities with bat 
SL-CoV ZC45, but less than 76.8% nt and 76.0% aa 
identities with the other members in the subgenus 
Sarbecovirus. In addition, ORF8 of the novel CoVs 
showed 88.5% nt and 94.2% aa identity with bat SL-CoV 
ZCA4S5, respectively, and less than 67.8% nt and 58.6% aa 
identity, respectively, with other members of Sarbecovirus. 
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Figure 2: Characteristics of viral genes. The schematic diagram of the novel coronavirus (CoV)’s genome (A). Phylogenic analysis of viral whole genome (B), spike (C), nucleocapsid (D) and 
RNA-dependent RNA polymerase (E) genes. The novel CoV identified and closely related viruses are in red in the phylogenetic trees. Other out-group viruses from public database are shown 
in blue. Evolutionary distances were calculated with the maximum likelihood method. Amino acid sequence of the putative receptor binding domain (RBD) of IPBCAMS-WH-01/2019 
compared with strains of severe acute respiratory syndrome coronavirus (SARS-CoV) and bat SARS-like CoV (F). 
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Table 2: Location and size of the putative proteins of the representative strain, IPBCAMS-WH-01/2019. 
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Similarity, % (nt/aa) 





Gene region Location“ Size, aa 
Whole genome NA NA 
S'UTR 1-265 NA 
nsp1 266-805 180 
nsp2 806-2719 638 
nsp3 (PLpro) 2720-8554 1945 
nsp4 8555-10,054 500 
nspS (3CLpro) 10,055-10,972 306 
nsp6 10,973-11,842 290 
nsp7 11,843-12,091 83 
nsp8 12,092-12,685 198 
nsp9 12,686-13,024 113 
nsp10 13,025-13,441 139 
nsp12 (RdRp) 13,442-13,468, 932 
13,468-16,236 

nsp13 (Hel) 16,237-18,039 601 
nsp14 (ExoN) 18,040-19,620 527 
nsp15 19,621-20,658 346 
nsp16 20,659-21,552 298 
S 21,563-25,384 1274 
3 25,393-26,220 276 
E 26,245-26,472 76 
M 26,523-27,191 223 
6 27,202-27,387 62 
7a 27,394-27,759 122 
7b 27,756-27,887 43 
8 27,894-28,259 122 
N 28,274-29,533 420 
9b 28,284-28,577 97 
3‘UTR 29,534-29,870 NA 


GC content, % SARS-Cov bat-SL-ZC45 
38.02 79.0 87.7 
44.53 88.7 93.5 
48.52 82.2/84.4 93.1/95.5 
39.86 68.3/68.3 92.2/95.2 
35.94 72.7175.8 91/93.8 
36.47 74.8/79.8 90.2/96.6 
38.13 83.6/96.0 92.1/99.0 
36.32 79187.2. 91.7/97.9 
38.15 79.9/98.7 82.7/100 
38.38 86.1/97.4 87.2/97.4 
40.12 84.9/97.3 86.7/97.3 
42.44 88.2/97.1 90.8/97.1 
37.28 88.5/96.3 86.5/95.9 
38.16 88.4/99.8 87.9/99.3 
38.27 83.3/95.0 82.9/94.4 
34.01 82/88.7 82.7/89.0 
35.46 85.1/93.2 89.5/97.9 
37.31 72.3175.5 75.9/80.4 
39.49 75.3/72" 87.8/90.9 
38.16 93.5/94.7 98.6/100 
42.6 85.2/90.5 93.4/98.6 
27.96 73.9/66.6 95.1/93.4 
38.25 82.1/85.2 88.7/87.6* 
31.06 82.0/79.5 NA 
35.79 NA 88.5/94.2 
47.2 88.1/90.5 91.1/94.2 
49.66 88.5/72.4 89.1/73.1 
40.36 95.2 93.2 





“Representative strain, PBCAMS-WH-01/2019, no. EPI_ISL_402123. Compared with 3a of SARS-CoV. * Compared with 7a of bat-SL-ZC45. aa: 
Amino acids; nt: Nucleotide; SARS-CoV: Severe acute respiratory syndrome coronavirus; SL: SARS-like; UTR: Untranslated regions; nsp: Non-structural 


protein; NA: Not applicable. 


These findings further indicated that the novel CoVs are of 
bat origin. 


The RBD in the CoV S protein determines the host 
range.''?! The RBD aa sequences of the novel CoV showed 
several distinct features, including higher aa identities with 
those of SARS-CoV (73.8-74.8%) and human angioten- 
sin-converting enzyme 2 (hACE2)-using SL-CoVs (76.4- 
76.9%) than those of SL-CoVs incapable of using hACE2 
(61.5-64.1%). The novel CoV does not possess the 
deletions commonly found in the RBD of SL-CoVs 
incapable of using hACE2 as a receptor [Figure 2F].'7°! 
In addition, the five critical aa residues interacting with 
hACE2 in SARS-CoV RBD (Y442, L472, N479, D480, 
T487) differ from the corresponding residues in the novel 
CoVs (L, F, Q, S, N), although these residues possess 
similar polarity.?°?"! These results suggested that the 
novel CoVs might still use hACE2 as the receptor. 


Viral culture 


CPE were observed in 30% of Vero cells inoculated with 
the new CoV after two passages [Figure 3A]. The cells 
showed a round, refractive, and syncytium appearance. 


The Vero cells with CPE were further examined using 
negative-staining electron microscopy, demonstrating 
characteristic CoV particles with surface projections 
[Figure 3B]. Immunofluorescent assays of the culture of 
Vero cells showing CPE with the convalescent serum 
from patients showed green signals in the cytoplasm, with 
no signals detected in wells containing control serum, 
indicating the presence of viral particles in the cells 
[Figure 3C]. 


Clinical features and outcomes of the patients 


The clinical features and laboratory test results of the five 
patients are summarized in Table 2. Fever, cough, and 
dyspnea were the most common symptoms. The white 
blood cell counts varied among these patients, but the 
lymphocyte counts were generally low. The alanine 
aminotransferase and serum creatine levels were normal 
or only slightly increased. Bilateral ground-glass opacities 
and consolidation were observed on chest radiography 
from two representative patients, Patient 2 based on aortic 
arch scan [Figure 4A] and pulmonary vein scan [Figure 4B] 
on day 10 after symptoms onset and Patient 5 taken on day 
12 [Figure 4C] and 13 [Figure 4D] after symptoms onset. 
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Figure 3: Viral isolation and identification. The untreated control (left) and cytopathic effect (right) on Vero cells visualized with an inverted microscope (A, original magnification x 20). Viral 
particles are negative stained with 1% solution of phosphotungstic acid and observed by using electron microscopy analysis (B, scale bar 200 nm). Immunofluorescent assays with 
convalescent serum and anti-human immunoglobulin G conjugated fluorescein isothiocyanate show signals (green) in cytoplasma in Vero cells. The Nuclei and cytoplasma were 
counterstained with 4’,6’-diamidino-2-phenylindole (DAPI, blue) and Evans blue (red), respectively (C, scale bar 20 wm). 


Several complications were observed in these patients. 
Four of the five patients (except for Patient 3) developed 
acute respiratory distress syndrome requiring oxygen 
therapy, and two patients were given extracorporeal 
membrane oxygenation. Two patients (Patients 1 and 5) 
experienced secondary infections, and Patient 5 later 
developed septic shock as well as acute kidney injury, and 
ultimately died of multi-organ failure. Patient 3 was 
discharged on January 8, 2020 (day 17 after symptoms 
onset). The other three patients were still hospitalized at 
the time of manuscript preparation. The treatments for 
these patients were shown in Table 1. 


Discussion 


In this study, we identified a previously unknown CoV 
from patients suffering from severe pneumonia. The 
whole-genome sequences of the viruses were obtained 
by a next-generation sequencing approach from all five 
patients, demonstrating overwhelmingly dominant viral 


reads in the BAL samples. Among the five novel CoV 
genome sequences, the nt identities reached up to 99.8% to 
99.9%. The viruses successfully isolated from the patients 
could also be effectively recognized by serum from a 
convalescent patient. These findings primarily indicate that 
the novel CoV is associated with the pneumonia that 
developed in these patients. However, it remains to be 
determined whether this novel CoV is capable of causing 
similar diseases in experimental animals. 


Sequence homology analysis of the viral genome showed 
that the CoV identified in this study is distinct from any of 
the known human CoVs, including SARS-CoV and 
MERS-CoV. The most closely related known viruses are 
two bat SL-CoVs (GenBank accession nos. MG772933, 
772934) identified in 2005 in Zhoushan, Zhejiang, China, 
which is geographically distant from Wuhan'**!; however, 
the nt identities among the viruses are only 85.7% to 
86.8%. Phylogenetic analysis showed that this virus forms 
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Figure 4: Computed tomographic chest radiographs of Patient 2, obtained on day 10 from illness onset at aortic arch (A) and pulmonary vein (B) scan demonstrating bilateral ground-glass 
Opacity and consolidation, and Patient 5 on day 12 (C) and 13 (D) after illness onset demonstrating white lungs. 


a single clade. Collectively, these data indicate that this 
CoV should be considered a new species. The outbreak of 
SARS in 2003 largely boosted awareness of threats caused 
by emerging CoVs. Consequently, great efforts have been 
made to monitor novel emerging CoVs and to trace their 
origins so as to establish a risk assessment and alert system 
for preventing potential epidemics in the human popula- 
tion. Clarification of the coronavirome in animals, 
particularly in bats as a key reservoir of a wide range of 
CoVs, should be a priority for any task force.?3-24) 


A few striking features of these novel CoVs indicated that 
they are of bat-origin. First, the genome sequences of the 
novel CoVs show high similarity with that of bat SL-CoV 
ZC45. Second, the phylogenetic analysis indicated that 
these viruses are evolutionarily close to bat SL-CoVs 
ZXC21 and ZC4S. Third, all of these novel CoVs contain 
ORF3 and intact ORF8 gene regions, which are 
characteristic features of bat-origin CoVs.!'77°! Moreover, 
the aa sequences of the N-terminal domains (NTDs) of the 
novel CoVs were very similar to those of ZC45 and 
ZXC21, whereas the RBD of the novel CoV showed higher 
aa sequence identity to that of SARS-CoV than to those of 


ZC45 and ZXC21, suggesting that a recombination 
event might have occurred at the region between the 
NTD and RBD of the S gene, facilitating the interspecies 
transmission. 


Owing to the lack of epidemic information at present, 
the transmission modes of the novel CoV remain obscure. 
It is notable that three of the five patients had a history of 
recent exposure to a seafood market in Wuhan. However, 
the origin of infection is unknown at the time of 
manuscript preparation. It is assumed that the zoonotic 
CoV jumped to humans through an intermediate host; 
for example, camel is suspected as the intermediate 
host of MERS-CoV, whereas the palm civet may have 
contributed to the interspecies transmission of SARS-CoV 
to humans.°-7°! Bat CoVs may evolve to adapt to 
using humans as a host during their circulation in a 
mammalian host, thereby enabling them to effectively 
infect humans,24 However, two of our patients did not 
have a history of exposure to the seafood market. 
Therefore, further investigation will be needed to deter- 
mine the potential of multiple infection sources responsible 
for this uncommon outbreak. 
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One of the most striking and concerning features of this 
virus is its ability to cause severe respiratory syndrome. The 
disease progressed rapidly with a major presentation of 
lower respiratory pathology. Notably, no obvious upper 
respiratory tract symptoms such as a sore throat and 
rhinorrhea were present in these patients. Therefore, 
further exploration is needed on the distribution of the 
viral receptor in the organs to potentially account for 
pathogenesis development. In addition, the possibility of 
unrecognized mild infections or subclinical infections 
should be clarified, as identification of such infections is 
critical to control spread of the disease. Development of 
serological assays would be largely beneficial to detect such 
types of infection at the population level. 


In conclusion, we identified a novel bat-borne CoV 
associated with a severe and fatal respiratory disease in 
humans. The emergence of this virus poses a potential 
threat to public health. Therefore, clarification of the 
source and transmission mode of these infections is 
urgently needed to prevent a potential epidemic. 
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